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Abstract. Recent results from the HERA experiment HI on the longitudinal 
stucture function Fl of the proton are presented. They include proton structure 
function analyses with particular emphasis on those kinematic regions which are 
^ ' sensitive to F^. All results can be consistently described within the framework of 

perturbative QCD. 
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^ ! 1. Introduction 



Deep inelastic lepton nucleon scattering (DIS) experiments have played a key role in the 
understanding of hadronic matter. They have revealed the structure of hadrons being 



^ ■ made out of constituents and the interaction between the latter. 

In fixed target DIS experiments scaling violations have been observed, i.e. the 
variation at fixed values of Bjorken-x of the structure functions with Q^, the squared 



four-momentum transfer between lepton and nucleon. These scaling violations are well 
described by perturbative Quantum Chromodynamics (pQCD). The evolution of the 
proton structure function F2{x,Q'^) is related to the gluon momentum distribution in 
the proton, xg{x,Q^), and the strong interaction coupling constant, as- Both, as and 
xg{x,Q'^), can be determined with precision ep DIS cross section data. 

HERA offers the possibility to study the structure of the proton over a wide 
kinematic range of x and Q^, several orders of magnitude larger than in the earlier 
fixed target experiments. Already first measurements of F2 at HERA, for x ~ 10~^ 
and ~ 20 GeV^, revealed a steep rise of -^2(0;, Q^) towards low x for fixed Q"^ IH, @]- 
Later measurements of F2 by H1[0, ||, Q and ZEUS0, 0] in a wide kinematic range, 
for values from 1 to 30000 GeV^ and for 3 ■ 10^^ < x < 0.65, have confirmed this 
trend. The strong scaling violations observed at low x are attributed to the high gluon 
density in the proton. The validity of the DGLAP evolution equation |p, which neglects 
higher-order terms proportional to -Inll/x), is questionable in the low x range 

and therefore has to be tested against the data. At sufficiently low x, non-linear gluon 
interaction effects have been considered in order to damp the rise of the cross section in 



accordance with unitarity requirements ||Tl . 



Determination of the longitudinal structure function Fl at HERA 



2 



The measurement of the longitudinal structure function Fl{x,Q'^) is of great 
theoretical importance, since it may allow distinguishing between different models 
describing the QCD evolution at low x. In fact, the structure function measurements 
at HERA remain incomplete until the longitudinal structure function Fl is actually 
measured. The direct, classical Fl measurement by substantial variation of the beam 
energies I T^] is yet to come at HERA. However, while waiting for this to happen in 



the HERA II period, the HI collaboration have extracted experimental results on Fl 
by measuring the cross section in a kinematic region where the Fl contribution is 
substantial (the high y region, see below) and comparing these measurements with 
the extrapolation of the F2 measurements. These results are described in this report. 



2. Neutral current cross section and the proton structure functions 

The cross section of the inclusive electron-proton DIS neutral current process (Figure 1), 
ep —> eX, depends on the three independent variables x and and s, the center of 
mass energy squared. and x are defined in terms of the four momenta of the incident 




Figure 1. Feynman diagram of neutral current DIS. 



and the scattered electron, k and k', and of the incident proton, P as 
Another important variable is the inelasticity 

which is related to x, and s hj = sxy Q. 

In the one-photon exchange approximation {Z^ exchange neglected: ^ '^lo) 
the cross section of the ep — * eX process can be expressed as 

da ~ L^,Wr (3) 

Here L^j^ denotes the leptonic tensor, describing the interaction between the electron 
and the virtual photon. In QED the lepton-photon vertex is well defined and therefore 

* Here and in the following formulae the proton and electron masses are neglected. 
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L^i, can be calculated exactly. The hadronic tensor W^^ corresponds to the photon- 
proton interaction and is not calculable. However, using Lorenz invariance and current 
conservation the hadronic tensor can be reduced to only two functions related to the 
structure of the proton. The double differential DIS cross section can then be written 

as 

One can consider the ep scattering process as the interaction of a flux of virtual 
photons [|l^ with the proton. The DIS cross section can then be written as 

^ = rfe)(.. + e(„..). (5) 

where r(?/) = Y^a/ {2t:Q'^x) stands for the flux factor, e(y) = 2(1 — y)/Y^ defines the 
photon polarization and 0"^, ai correspond to the cross sections of the interaction of 
transverse and longitudinally polarized photons, respectively. These cross sections are 
related to the structure functions as 

A.^a (6) 

The longitudinal structure function Fl is directly proportional to o"l(x, Q"^). Due to the 
positivity of the cross sections, the two structure functions F2 and Fl obey the relation 

< < F2. (7) 

The "reduced" cross section is defined as 

r_, 

Since the contribution of the longitudinal structure function F^ to the cross section can 
be sizeable only at large values of in a large kinematic range the relation Or ~ F2 
holds to a very good approximation. 

In the Quark Parton Model (QPM) [l^ the structure function F2 can be expressed 
as a sum of the quark-antiquark momentum distributions xqi{x) weighted with the 
square of the quark electric charges 

F2{x) = J2 e- a;(gi(a;) + g,(a;)). (9) 

i 

For spin 1/2 partons QPM also predicts ctl = 0, which leads to the so-called Callan- 



ar = F2ix,Q^)-^FLix,Q^). (8) 



Gross relation |]15|: Fl{x) = 0. 

The naive QPM has to be modified in QCD as quarks interact through gluons, 
and can radiate gluons. Radiated gluons in turn can split into quark-antiquark pairs 
("sea quarks") or gluons. The gluon radiation results in a transverse momentum 
component of the quarks. Consequently, quarks can also couple to longitudinally 
polarized photons and the Callan-Gross relation is no longer satisfied exactly. Thus, in 
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QCD the longitudinal structure function is non-zero. Due to its origin, F^ is directly 
dependent on the gluon distribution in the proton and therefore the measurement of F^ 
provides a sensitive test of perturbative QCD. 



3. Experimental procedure 
3.1. Kinematic reconstruction 

At HERA 27.5 GeV positrons collided with 820 GeV protons. For the last two years 
the proton beam energy was increased to 920 GeV, and consequently the center of mass 
energy from 300 to 319 GeV. 

The kinematic variables y (or x) and can be reconstructed either from the energy 
and the polar angle of the scattered positron, or from the energy and the polar angle 
of the hadronic final state, or using a combination of these energies and polar angles. 
The overconstrained kinematics allows cross checks, as well as the choice of the method 
with the best experimental precision for a given kinematic range. 

In the "electron method" the event kinematics is reconstructed using the energy 
E'^ and the polar angle Oe of the scattered positron 

Ql = ye = l- ^s^n\eJ2). (10) 

I - ye Ee 

While the electron method is accurate at large values of y, corresponding to low energies 
E'^, the resolution rapidly degrades for y^ where E'^ approaches the positron beam 
energy Ee. The kinematic variables at low y are reconstructed using also information 



from the hadronic final state. The inelasticity y can also be determined as [16 



'^' = 2Ee = 2^;' (^^^ 

where Ei and pz,i are energy and longitudinal momentum components of the particle i 
in the hadronic final state. In the "S method" []T^ and y are reconstructed as 



2 ^ E'^sin'^e^ ^ yh ^ .^2) 

The S method has good resolution also at low y, and is insensitive to the initial state 
radiation. The cross section is extracted using a combination of the electron method 
and the S method. 



3.2. HI detector 

The HI detector [l^ combines tracking in a solenoidal magnetic field of 1.15 T 
with nearly hermetic calorimetry to investigate high energy ep interactions at HERA. 
For polar angles 6e > 153° the scattered positron energy E'^ is measured in the 
backward electromagnetic lead-fibre scintillator calorimeter (SPACAL) [^] which has an 
integrated timing function to veto proton beam induced background interactions. Here 
9 is defined with respect to the proton beam direction. Identification of the scattered 
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positron is improved and the polar angle measured with a backward drift chamber 
(BDC), situated in front of the SPACAL, and with the new backward silicon detector 
(BST) The BST consists of four detector planes, arranged perpendicular to the 

beam axis and equipped with 16 wedge shaped, double metal silicon strip detectors. 
The BST measures the polar angle of the crossing tracks with an internal resolution of 
about 0.2 mrad for 9 between 172° and 177°. For polar angle A° < 9 < 154° the scatterd 
positron is measured in the LAr calorimeter. 

The hadronic final state is reconstructed using the Liquid Argon (LAr) calorimeter, 
tracking detectors and the SPACAL. The interaction vertex is determined with the 
central drift chambers, the jet chamber CJC and two layers of z drift chambers, mostly 
using the hadronic final state particles. This vertex determination is complemented 
by the inner proportional chamber CIP , for 167° < 9^ < 171°, and by BST, for 
171° <9e< 176.5°. 

The luminosity is determined with a precision of 1.5% using the small- angle 



bremsstrahlung process ep ep7 The final state photon and positron, scattered 
at very low Q^, can be detected in calorimeters ("photon and electron taggers") which 
are situated close to the beam pipe at distances of 103 m and 33 m from the interaction 
point in the positron beam direction. 



3.3. Event selection and simulation 

The event selection criteria are slightly different in the low and in the high ranges. 
In the low range the scattered positron is identified in the SPACAL. In the high 
range the positron is seen in the LAr calorimeter as that cluster of maximum transverse 
energy, for which the requirements on the cluster shape are satisfied. The difference is 
then mainly in the cluster definition and in the track validation. In the LAr cluster-track 
validation we use CJC tracks for 9 > 35°, in the SPACAL cluster case we use CJC, BDC 
and BST tracks, where appropriate. 

Longitudinal momentum conservation in neutral current DIS events gives the 
constraint that E — Pz, summed over the final state particles, is about equal to 2Ee. In 
events with initial state radiation the radiative photon may carry a significant fraction 
of the E — Pz sum. Such events are thus suppressed by requiring E — pz > 35 GeV. 

Distributions of the selected events were compared with simulations of deep inelastic 
scattering, photoproduction and Compton scattering. The simulated Monte Carlo 
events were submitted to the same reconstruction and analysis chain as the real data. 
The track detection efficiencies were determined from the data utilizing the redundancy 
of the central and backward tracking detector system. Agreement at the few per cent 
level was reached between experiment and simulation for the large set of technical and 
physics distributions studied. 



4. The longitudinal structure function Fl(x,Q'^) 
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Photoproduction background estimate 

For the measurement of the longitudinal structure function it is essential to reach the 
highest possible values of y. This requires the identification of scattered positrons with 
energies down to only a few GeV as well as the efficient rejection of photoproduction 
background events in which low energy deposits in the SPACAL can mimic the signature 
of a deep inelastically scattered positron. 

The main part of the background is due to photons from tt° — 77 decays. In 
the low region a sizeable fraction of this background can be removed by requiring 
a track signal in the BST. The remaining background is due to photon conversion and 
showering in the passive material in front of the detectors, possible overlap of n" decays 
with charged tracks, and misidentified charged pions. It is subtracted bin by bin using 
the PHOJET p2| MC simulation. The photoproduction background can be estimated 
experimentally using tagged events in which the scattered positron is detected in the 
electron tagger. 

At above 10 GeV^, for y < 0.75 the photoproduction background is subtracted 
using the PHOJET simulation. For y values above 0.75 experimental information is used 
by employing the charge assignment of central tracks associated with SPACAL energy 
clusters. This allows the energy range to be extended down to E'^ =3 GeV, corresponding 
to y < 0.89. The candidates with negative charge are taken to represent the background 
in the positron data sample. This statistical subtraction procedure requires the study 
of any process which may cause a charge asymmetry. Suchs asymmetry is e.g. due to 
the antiproton interaction cross section exceeding that for proton interactions at low 



energies p3|. Annihilation leads to larger energy deposits in SPACAL than proton 
interactions which introduces an asymmetry for low energies above a given threshold. 
This charge asymmetry can be measured using tagged photoproduction events which 
fulfill the DIS event selection criteria and is taken into account in the measurement of 
the positron DIS cross section at high y. 

As already stressed above, the detailed understanding of the data down to 3 GeV 
scattered positron energy is essential for the measurement of the DIS cross sections at 
high y, enabling the determination of Fl. Control plots illustrating the cross section 
measurements at high y are shown in Figures 2 and 3. 



4.2. Cross section at large y 

The cross section is determined by converting the measured number of events, after 
background subtraction, in bin averaged cross sections using acceptance calculations 
from the Monte Carlo sample. The precision for this measurement is dominated by 
systematic uncertainties of typically 3%, extending to about 7% at the edges of the 
covered y range. 

The measured reduced cross section cXj. is shown in Figures 4 and 5 for two e^p 
data samples, taken in 1996/97 and 1999 with 820 GeV and 920 GeV proton beam 
energy, respectively. The cross section rises towards low x, turning over at x values 




Figure 2. Distributions illustrating the cross-section measurement at high y 
(0.46 < y < 0.82) and 2 < < 5 GeV^ for events in the BST acceptance 

range. DIS event distributions of a) the polar angle and b) the SPACAL energy of 
the scattered positron, c) SPACAL energy distribution for tagged photoproduction 
events fulfilling the DIS event selection criteria, apart from the E ~ requirement. 
Solid points: HI data; shaded histograms: simulation of photoproduction events; open 
histograms: added distributions of simulated DIS and photoproduction events. 
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Figure 3. Distributions illustrating the cross-section measurement at high y 
(0.46 < y < 0.89) and 10 < < 35 GeV^. a) Polar angle and b) SPACAL 

energy distributions before subtraction of the photoproduction background using the 
charge measurement by the CJC. Solid points: data with positive charge assignment. 
Shaded histogram: data with negative charge assignment. Open histogram: sum of 
data with negative charge assignment and DIS event simulation, normalised to the 
integrated luminosity of the data, c) Spectrum of energy measured in the electron 
tagger for DIS candidate events with a linked track of either positive charge (solid 
points) or negative charge (histogram). 



corresponding to y ~0.6. This behaviour, for > 3.5 GeV^, is well described by the 
NLO QCD fit 1^ to the 820 GeV data. Note that this fit used HI data with a minimum 
Qmin ~ 3-5 GeV^. The extrapolation of the fit to values lower than 2.5 GeV^ falls 
below the measured data. 
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4-3. Determination of Fi 

The extraction of the longitudinal structure function Fl is based on the reduced cross 
section (8), which depends on F2{x,Q'^) and Fl^XjQ"^). The effect of Fl is enhanced 
with y"^, and the reduced cross section tends to F2 — F/, as ?/ — 1. For the extraction 
of the longitudinal structure function from the inclusive cross section the data at large 
y are very important. An important advantage of HERA, compared to fixed target 
DIS lepton-nucleon experiments, is the wide range of y values covered. This allows the 
behaviour of F2 at low y to be determined reliably and to be extrapolated into the region 
of high y. 

Two methods are used by HI to perform the extraction of the longitudinal structure 
function: 

For > 10 GeV^ the "extrapolation method", introduced in ||25|, is used. A NLO 



DGLAP QCD fit is used to extrapolate F2 into the high y region. This fit uses only HI 
data in the restricted kinematic range y < 0.35 and > 3.5 GeV^. In Figure 4 the fit 
is compared with the measured cross section for bins above 10 GeV^, accessing the 
high y region. The difference between the measured ar and the extrapolated F2 is used 




Figure 4. Measurement of the reduced DIS cross section (closed points). Triangles 
represent data from the BCDMS muon-proton scattering experiment. The curves 
represent a NLO QCD fit to the HI data alone, using data with y < 0.35 and > 3.5 
GeV^. The dashed curves show the F2 structure function as determined with this fit. 
The error bands represent the experimental and model uncertainty of the QCD fit. 
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Figure 5. Measurement of the reduced DIS cross section (closed points) from the 
e^p data taken in 1999 at 920 GeV proton energy |Q. The curves represent a NLO 
QCD fit to the HI 1996/97 data, which was hmited to > 3.5 QeV^. The dashed 
curves show the extrapolation of this fit below Q^„j„. For each bin the lower curves 
(turning at low x) represent the reduced cross section and the upper curves represent 
the structure function which rises towards low x. For large x, corresponding to 

low y, CTr — -^2- 



to determined Fl{x,Q'^). Systematic errors, which are common to the lower y and the 



larger y region, are considered in the fit as described in |26l . 

At low < 10 GeV^, the behaviour of F2 as a function of Iny is used in a new 
extraction method p. This so called " derivative method" is based on the cross section 
derivative {dar / d\n.y)Q2. 

The derivative of the reduced cross section, taken at fixed Q^, is given by 

dar.\ f dF,\ ^ 2-y OFl y' 



d\ny J Q2 \d\ny J Q2 Yf d\n.y Y+ 

For y — > 1 the cross section derivative tends to the limit {dF2/dlny)Q2—2-FL, neglecting 
the contribution from the derivative of Fl. At largest y the Fl contribution dominates 
the derivative of the reduced cross section cr^- This is in contrast to the influence of 
Fl on which is dominated by the contribution of F2 for all y. A further advantage 
of the derivative method is that it can be applied down to very low ~ 1 GeV^ 
where a QCD description of F2{x,Q'^) is complicated due to higher order and possible 
non-perturbative corrections. The cross section derivatives are shown in Figure 6 as 
function of y. This measurement of dar/dlny was used to determine the longitudinal 
structure function. The assumption has been made that the derivative dF2/dlny is 
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Figure 6. Measurement of the derivative {dar / d\iiy)Q2 . The curves represent the 
QCD fit result to the HI 1996/97 data for y < 0.35 and > 3 5 QeV^ calculated 
with different assumptions on Fl- The inner error bars represent the statistical errors 
and the total error bars the statistical and systematic errors, added in quadrature. 



a linear function of Iny up to large y. In each bin straight line fits were made 
to the derivative data for y < 0.3. The line fits describe the data very well and the 
extrapolation of the straight line was taken to represent the contribution of F2 at high 
y. The extrapolations were compared with the values obtained from QCD and good 
agreement was found. The uncertainties of the straight line extrapolation were included 
into the systematic errors of the measurement. The small contribution of dFi/dhiy to 
the derivative was corrected for by using NLO QCD and the size of the full correction 
was added to the overall error of the measured F^ . 



4-4- Results 

The longitudinal structure function Fl{x,Q'^) determined by the derivative method is 
shown in Figure 7. The data extend into the lower x range and are consistent with 
previous HI results. The curves represent NLO QCD calculations of Fl based on a fit 
to the 1996/97 data for y < 0.35, i.e. in a region where F2 can be measured practically 
independently of Fl- 

In Figure 8 the determination of Q^) at high is shown for e^p and e~p 
data sets. The extrapolation method is used. Both data sets are consistent and in 
agreement with the QCD fit. The possible extreme values for Fl [Fl = and Fl = F2) 
are shown for comparison. 

An overview of all current HI data on Q^), from Q"^ = 2.2 GeV^ to 700 GeV^ is 
given in Figure 9. The data extend the knowledge of the longitudinal structure function 



Determination of the longitudinal structure function Fl at HERA 



11 



;0.75 

0,5 

0.25 


1 

0,75 
0,5 
0,25 




Q'=2,2 GeV' 



Q'=4,2 GeV' 



i HI 99 prel, 
• HI 96-97 
-QCDfit(HI) 




Q'=7,5GeV' : Q'=12GeV' 





-5 -4 -5 -4 

10 10 10 





n 







10 



Figure 7. The longitudinal structure function Fl{x,Q^) as a function of x in bins 
of Q^, determined with the derivative method from the e^p data taken in 1996/97 
(points) and in 1999 (triangles). The inner error bars represent the statistical error 
and the total error bars the statistical and systematic errors, added in quadrature. 
The curves represent the QCD fit result to the HI 1996/97 data for y < 0.35. 



into the region of low x, much beyond the region of fixed target lepton-proton scattering 
experiments. The increase of Fl{x, Q^) towards low x is consistent with the NLO QCD 
calculation, reflecting the rise of the gluon momentum distribution in this region. The 
values of Fl{x,Q'^) are thus severely constrained by the present data, unless there are 
deviations from the assumed extrapolation of F2 into the region of large y corresponding 
to the smallest x. A measurement of the x dependence of Fl{x,Q'^), independent of 
assumptions about the behaviour of F2, can be performed with a variation of the proton 
beam energy at HERA 



A NLO QCD fit has also been performed by the ZEUS collaboration [27]. Within 
this fit Fl was calculated, and found to be consistent with the HI data. 

Several theoretical models were compared to the HI data. Predictions based on 
the saturation model |^ are shown in Figure 10 together with the HI Fi{x,Q^) data 
points. At large x and values data are well described by the model, while at low 
X and low higher twist contributions are clearly needed to describe the Fl{x,Q^) 
distribution. 

In Figure 11 the HI data are compared to the Donnachie-Dosch model 



predictions. This model is based on a dipole picture in which large dipoles couple 
to a soft pomeron and small dipoles couple to a hard pomeron. The parameters are 
fixed by proton-proton scattering data and by the measured proton structure function 
F2{x,Q'^). This simple approach is able to describe the Fl data well. One can see 
that with increasing the importance of the hard pomeron is increasing, and it is 
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Figure 8. The longitudinal structure function Fl as a function of at fixed y =0.75 
for the e^p and e~p data. The inner error bars represent the statistical error and the 
total error bars the statistical and systematic errors, added in quadrature. The curves 
represent the QCD fit result calculated with different assumptions on F^. The shaded 
band shows the expectation for Fl, and its uncertainty, from the QCD fit. 



dominating at the high values. 
5. Summary 

New results on the extraction of the longitudinal structure function Fl{x,Q'^) are 
presented for a wide range from 2.2 GeV^ to 700 GeV^ and for much lower x 
values (3 • 10^^ < x), compared to the fixed target experiments. The extraction results 
are possible due to the improved detectors in the HI backward region, used for the 
identification and measurement of the scattered positrons at low and high y. 

Two different methods are used to determine F^. The partial derivative of the 
reduced cross section, {dar/d\iay)Q2, is used to extract Fl at low < 10 GeV^. At 
> 10 GeV^ the difference between the measured reduced cross section 0"^ and F2, 
calculated from an extrapolation of a NLO QCD fit to low y data, is used for the Fl 
determination. Thus the longitudinal structure function Fl at low x is determined more 
precisely than hitherto and in a larger range. 

The observed rise of FlIxjQ"^) towards low x is consistent with the NLO QCD 
calculations and reflects the rise of the gluon momentum distributions in this kinematical 
region. 

The data are also well described by calculations based on a colour dipole model 
and on a combination of the colour dipole with a two pomeron model. 
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Figure 9. The longitudinal structure function Fl{x,Q^) as a function of x in bins 
of Q^, obtained by HI and by charged lepton-nucleon fixed target experiments. The 
inner error bars represent the statistical error and the total error bars the statistical and 
systematic errors, added in quadrature. The error bands are due to the experimental 
(inner) and model (outer) uncertainty of Fl using the NLO QCD fit to the HI 1996/97 
data for y < 0.35 and > 3.5 GeV^. 
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Figure 10. The longitudinal structure function Fl{x, Q^) as a function of x in bins 
of Q^, obtained by HI and by charged lepton-nucleon fixed target experiments. The 
inner error bars represent the statistical error and the total error bars the statistical 
and systematic errors, added in quadrature. The dash-dotted and dotted lines are 
fits using the saturation model of Golec-Biernat and Wiisthoff , with higher twist 
contributions (exact) and without higher twists (LT). 
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Figure 11. The longitudinal structure function as function of x in bins of Q^i as 
measured by HI. Curves represent the Donnachie-Dosch model |p9| predictions. The 
solid line is the full result and the dashed line is the hard pomeron contribution. 
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